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Abstract 

From SUSY ladder operators in momentum space of a neutron in 
the magnetic field of a linear current, we construct 2x2 matrix oper- 
ators that together with the z-component of the angular momentum 
satisfy the su(2) Lie algebra. We use this fact to explain the degener- 
acy of the energy spectrum. 

PACS: 11.30.Pb; 03.65.Fd 



1 Introduction 

SUSY QM is the modern way to study solvable and perturbative systems as 
is extensively shown in Refs. [1-5] . Recently, it has been shown that for some 
solvable systems in more than one dimension, the SUSY operators can be 
used to explain the degeneracy of the spectrum [6-10]. 
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The existence of magnetically bound states of a neutron embedded in the 
magnetic field of a filamentary current has been demonstrated experimentally 
by Schmiedmayer [11] and Vestergaard et. al. [12]. In the recent years this 
problem became very important from the experimental point of view. It is 
considered the basis for the so-called neutron magnetic storage, for the ultra 
cold neutrons physics and for the generation of radio frequency radiation [T3] . 

Pron'ko and Stroganov [H] were the pioneers in studying this system 
analytically in momentum space. Even though they explained the degeneracy 
of the system, the origin of the symmetry generators is not given. Bliimel 
and Dietrich [151 US] studied the problem and found the energy spectrum and 
the eigenf unctions in configuration space. However, the works [14-16] do not 
show the supersymmetry of the problem. In [T7] Voronin solved the problem 
in momentum representation from SUSY QM. By expanding the spinor wave 
function in Fourier series, and after some transformations, Voronin obtained 
an exact supersymmetric form for each component of the expansion. He also 
found the energy spectrum and the eigenf unctions of the system. We note 
that the theoretical description developed in [12J, and more recently in [18], 
was based upon the concept of SUSY QM for multicomponent wave functions 
in configuration space. 

In this Letter we study the problem in momentum space. We construct 
its symmetry group generators from SUSY QM and explain its degeneracy 
of the energy spectrum. 

Using polar coordinates we obtain an exact supersymmetric pair for the 
complete radial components of the spinor wave function. The second-order 
differential equation for each radial component is obtained, which is identi- 
fied with the first Poschl- Teller equation. By using the SUSY QM approach 
[19,20], wwe get the ladder operators. Generalizing the supersymmetric op- 
erators to matrix operators of two variables, it is shown that these operators 
together with J z satisfy the so(3) or su{2) Lie algebra. Using angular mo- 
mentum theory, we find the energy spectrum, and show that these matrix 
operators allow us to find the set of states correspondng to a given energy 
level. 

The magnetic field of an infinite straight wire earring a current / located 
along the z-axis 

g _ fj.pl (-y,x,0) ( 
2tt (x 2 + y 2 )' { ' 

Thus, the Hamiltonian of a neutron of mass M = 1 and magnetic moment 
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fia interacting with B is 



W = £ + I *.B = £ + G=ff±^ (2) 
2 2 (x z + y z ) 

where G is a constant defined as 

G=^. (3) 

Z7T 

and we have set c = h = 1. Considering the translational symmetry along 
the z-axis, the two component wave function of the system can be written as 

^ {rl= i 2mkz/L (Mx,y)\ (4) 
v ; y/z \M x >y)J 

where k — 0, 1, 2, 3, ... As a consequence of the free motion along the z— axis, 
we get a two dimensional problem. 

Thus, the Schrodinger equation can be written as a system of two differ- 
ential equations 

Px r, Py ^i{ x ^v) +iG — lV J ) 2 {x,y) = Efafay), (5) 



x z + y z 

,2 -2 

'h.(x, v) — iG 

x z + y z 



Py -ij 2 (x,y) -iG X 2 + ^ y 2 i)x(x,y) = Eif) 2 (x,y), (6) 



where E = E — 



L 2 ■ 



It is easy to show that 



H, J z 



0, where J z = L z + s z is the ^-component 

of the total angular momentum, therefore, using polar coodinates (p, varphi), 
((!]) can be written as 

where j z = ±|, ±|, ... is the quantum number corresponding to the operator 

If we multiply Eqs. (J5J) and (EJ) by (x + iy) and (x — iy), respectively, we 
are able to write these equations in momentum space in the form 

d d \ 

+ i— aip[ z (p x ,p y ) = -G^i z (p x ,p y ), (8) 



dp x dp. 
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d d 
— i- 



()j> "Qp J vfy'iPxiPy) = G^{ z (p x ,p y ), (9) 

where the variables x and y have been substituted by their corresponding 
operators in momentum space. In Eqs.flBD and we have defined 



as *±£-E = £-E (10) 



and ipa(Px,P y ), a = 1, 2 are the Fourier transforms of the functions ip a (Px,Py) 
given by 

i /*oo /»oo 

Mp X iPy) = —l I M^y)e- i{xpx+yPv) dxdy (11) 



oo J — oo 



Using polar coordinates in momentum space (p, 8), we write the cartesian 
components of the momentum p as follows 

p x =p cos 9, (12) 

p y = psin9. (13) 
Thus, Eq.f llip can be written as 

-i /*oo p2n 

#"(P'*) = 2^y o I ^ z (p^),e-^ pco < d -^pdpd V (14) 
where ipl z (p,(p) are defined in Eq.([7j). Therefore, 

poo 

$«(p, 0) = r^) e iO*Ti)0 / pfi(p)J, T i(pp)dp = e^^Ftip), (15) 

where Jj zT i are the Bessel funtions of order j z =)= |. The result given by 
equation ( fl5i) permits to write the spinor wave function as 

^+^Ff(p)J' (16) 

which implies the conservation of the operator J z in momentum space, as 
was expected. Explicitly, we have 



Therefore, Eqs. flHD arid are transformed to 

(4- - — ) aF tiP) = -GFt(p), 
\dp p J 

If we perform the change of variables defined by 

q=p(-2E)- l '\ 

Zt=p 1/2 (aFn 
now, Eqs. (ITB1) and ( fT9l) can be expressed as 



respectively, with 



qj \ + 1 

dq q J f + 1 



xl/2 - )fT ' 



h(-2Ey/ 2 



Now, setting g = tan |, we introduce the angular variable /3, < 
Thus, equations (1221) and (1231) are transformed to 



where 

A 

W(0 = ^(tanf + cot§ 



If we consider the following definitions, 



A< = -- + W W , A = - + W W , (29) 

Eqs.f l25|) and ff26l) can be written as 

ftZ{' = e l ' 2 Z 3 2 % (30) 

AZ 3 2 ' = eV 2 Zi'. (31) 



We recognize that Eqs. (l25l) and (1261) are written in the so-called exact su- 
persymmetric form, being W((3) the superpotential [20]. Operators A 1 and 
A can be used to define the supersymmetric Hamiltonians Hj z .+ = AA* and 
H jz ._ = A^A. From Eq.flUD and ([HID we can see that they operate on Z{ z 
and Z 2 , respectively. Say 

i^ :+ Zf = iUtzf' = eZf', (32) 

H jz: -Z£ =tiAZt =eZ£. (33) 

This means that the radial components of the spinor (|53j) . Z{ z and Z^ z , are 
the eigenfunctions of the SUSY Hamiltonians Hj z:+ and Hj z: _, respectively. 
Notice that, from fl32l) and fl33l) . the spectrum of the Hamiltonians Hj z -.+ and 
Hj x :- is the same e = e*, which implies a broken supesymmetry |21j . 

It is mmediate to find the spectrum e because we have identified ( 1281) as 
the superpotential of the problem. By using Table I of Ref . [20] , we find that 
the superpotential fl28l) is a particular case of 

W = A tan ax — B cot ax, (34) 

which corresponds to the Poschl- Teller I potential, with spectrum 

= (A + B + 2na) 2 -(A + B) 2 (35) 

and the ground state 

^ = (sinax) B/a (cosax) A/a . (36) 
For our case, if we set A = —B = j z /2 and a = 1/2, then 

e ± = n 2 , n = 1,2,3,..., (37) 
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ipQ = (sinaa;) 3 '* (cos ax) Jz . (38) 

Because, j z = ±|,±|,..., the ground state (138|) is non-normalizable. It 
means that the state with n = in e~ is abscent, which implies = e + , 
as we have moted after equation (I3"3"|) . This confirms that SUSY is broken 
|21j . It is remarkable that we have used changes of variable analogous to 
those used by Voronin[17j and obtained similar results. Nevertheless, the 
difference between our and that of Voronin is that we have not expanded 
the wave function in space. This allowed us to find the SUSY pair (125]) and 
( T25|) for the complete radial wave funtion in momentum space. The main 
results we will obtain in this Letter are based on the following fact: since 
the supersymmetric Hamiltonians Hj z .+ and Hj z .- are differential operators 
of second order, each of them can be factorized in a different way to SUSY. 
Indeed, by the Infeld-Hull (IH) factorization method [19]. Explicitly, Eqs. 
fl32]) and (J33D are 

H^tM = (-§2 + V ^m) = (39) 

with 

(Jz ~S Z + \){jz ~S z -\) (J z + Sz + \)(Jz + S Z ~ |) 1 

■20 0/3 ' ^ ' 

sin | cos^ | 

where we have defined = Z{ z {{3) and = Z[ z {f3). According to [19], 
we find that the IH ladder operators corresponding to tje first Poschl- Teller 
potential (f4"U|) are 

B S ,M - 

where 

By straightforward calculation, 



= -^ + M«, (41) 

= ^ + M«, (42) 

/ 2 V 2 4 4/ 2 v ; 

it is easy to show that the IH ladder 
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operators factorize the Hamiltonians Hj x .± in the following way 



Hi 

Jz 



(44) 



H 



Eqs. (1441) and (1451) can be written in the matrix form as 



H jzr 



o \ (*m 



(45) 



(46) 




H 



Jz- 



{(3)B Jz+l , + {(3) 






B) z+ ^_{(3)B 3z+1 ,.{(3) 




*7M 



(47) 



We have found the IH one variable ladder operators Bj Zt ± and Bj z . ± which 
factorize the corresponding Hamiltonians Hj x .±. Notice that these operators 
depend on the quantum number j z . To consturct the symmetries of the 
system, we generalize the IH ladder operators and Hj z; ± changing j z by its 
corresponding operator in the variable 9. This is achieved as follows. From 
Eqs. (|32|) and f|33l) . and the definitons after Eq. (1401) . Hj z - + and Hj x .- act on 
the upper (<&+(/?)) and the lower ($"(/?)) radial components of the spinor 
wave function, respectively. According to Eq. (IT7j) . the quantum number j z 
must be changed by + | if ti is within an operator acting on $+(/3), and 
by — i-Kz — | if ti is within an operator acting on $~ (j3). Therefore, operators 



8 



defined in Eqs.( l4"Tj) and (f4"21) can be generalized to 



(_d_ 

1 dp + 2[\ l de~2j 



cot 







1 



tan 



(48) 



-(9 



_ A I 

"d(5 + 2 



. d 1\ /? 
-z — — I — cot — 
06 2 2 



P 



(49) 



and the Hamiltonians #,v± are generalized to 



H 4 



1 

<9/3 2 + 4 



5 



5 



1 



06 







+ -* ™ - ^ -* ™ + r esc" - 



50 



/3 



-/— + 1 + - 1 I -z'^r ± 1 - - ] sec^ - 
06 2 J I 06 2} 2 



(50) 



Notice that in Eqs. fl48|) and (j49jl we have added a phase factor because 
these operators act on the omplete spinor components e l<Jz± 2)*j z (^). From 
Eqs.f l48p . (1491) and (|50|) we define the 2x2 matrix operators 













(• 




|. ■»-) 





51/ 



(51) 



and 



H 




Also, we define the spinor wave function 



(52) 



(53) 
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The definitons given by Eqs. (l5"Tj) . 
matrix equations (jlBjl and (pT7|) as 



and (1531) allow us to generalize the 



Lj+Zj^M) 



H 



J z + 



(54) 



(55) 



The last two equations can be written in compact form as 



J ± J T Z jz (P,9) 



H — I J, ± 



1 



(56) 



These equations, together with 

J,Z ia {p t 9) = (4 + s z ) Z jz ((3, 9) = Jz Z Jz (f3, 9), 
lead us to the commutation relations 



Jz, J± 



Z j ,{p,0) = 2J,Z jm {p, 
Z Jz (P,9) = ±J ± Z 3z (P, 



(57) 



(58) 



(59) 



which is the Lie algebra of one of the groups £0(3) or SU(2). By using this 
angular momentum algebra, we find that 



j 2 J. 



J 2 J Z 







and 

j T j ± = J 2 -J Z (J Z ±D* 

Also, from Eqs. (156]) and (!60|) it follows that 

[J 2 ,H} = 0. 



(60) 
(61) 

(62) 



Thus, from the IH factorization method we have constructed the operators 
J±. These operators together with J z allowed us to obtain the Casimir opera- 
tor J 2 , which commutes with J z and with the Hamiltonian H. Therefore, H, 
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J 2 and J z are a compkete set set of commuting operators which have simul- 
taneus eigenfunctions. This fact permits to introduce the principal quantum 
number j which satisfies 

j 2 Z jjK (P,8)=j(j + l)Z j j K (/3,8) (63) 

J t Z jdx {p,0)=j t Z jd ,(p,e), (64) 
with — j < j z < j. Substitution of (IBTj) into Eq. (156]) leads us to 

2" 



J — J z [ J z ±1 



H 



J z ± 



1 



Z J)ja {p,0). (65) 



From each of these equations we find that the spectrum of the operator H is 
given by 



(66) 



which is in agreement with that reported in the literature [22]. Using Eqs. 
(EH), (EZD and (JSSJ), we find that the spectrum of H is 



G 2 



, 2 ' 



(67) 



2 (* + §)' 

From the commutation relations (|59p . it is immediate to show that 

j z j±Z j:jz = (j z ± 1) J±Z j>jz , (68) 

which implies that j±Zj j z is an eigenvector of J z with eigenvalue j z ±l. This 
means that J±Zjj z oc Zj j z ±i. In fact, by taking the same convection of the 
angular momentum theory, we get 



.LZ JJZ+1 (P,9) = y/j(j + i)-j,(j z + 1^.(13, 



(69) 
(70) 



Notice that Eq. (T27j) was obtained by using SUSY QM, whereas Eq. (|66|) was 
obtained using angular momentum theory. Regarding that j z = ±|, ±|, ±|, .. 
then, j = |, |, |, thus, Eq. fT27]) and fl66l) reduce to the same energy spec- 
trum. Since each of these equations does not depend on the quantum number 
j z , the energy spectrum is degenerated. The restriction — j < j z < j implies 
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that the degeneracy of the energy spectrum ej is 2j + 1. This number is not 
easy to find by other way to that given here. A very important result of our 
treatment is that, because j = |, |, |, the symmetry gropu of the problem 
is SU(2), the covering group of the rotation gropu 5*0(3) [23] . 

Eqs. (1691 and flTOl show that operators J± can be used to obtain the set 
of degenerate states belonging to the enery level €j. 

We note that the energy spectrum, Eq. (I6TI) . as well as its degeneracy 
are equal to thoseof the two-dimensional hydrogen atom [24,25]. Moreover, 
it has been shown that, as a limit case, the degenerate wave wunctions and 
the spectrum of the problem we have considered here, reduce those of the 
two-dimensional hydrogen atom [T5], whose symmetry group is SO (3) [25] . 

We have obtained the symmetries of the system in momentum space 
bacause it is a two-dimensional system. This procedure os very difficult to 
perform for systems in more than two dimensions [TO] . 

As a final remark, the procedure described here can be applied to the 
same problem in configuration space. This is the subject of a forthcoming 
report. 
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